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ABSTRACT 

This paper presents the algorithms created and implemented by the 
Goddard Space Flight Center’s (GSFC’s) Solar Maximum Mission (SMM) 
attitude operations team to support large-angle spacecraft pointing at sci- 
entific objectives. The mission objective of the post-repair SMM satellite 
was to study solar phenomena. However, because the scientific instru- 
ments, such as the Coronagraph/Polarimeter (CP) and the Hard X-ray 
Burst Spectrometer (HXRBS), were able to view objects other than the 
Sun, attitude operations support for attitude pointing at large angles from 
the nominal solar-pointing attitudes has been required. Subsequently, atti- 
tude support for SMM has been provided for scientific objectives such as 
Comet Halley, Supernova 1987A, Cygnus X-l, and the Crab Nebula. In 
addition, the analysis has been extended to include the reverse problem, 
computing the right ascension and declination of a body given the off-Sun 
angles. This analysis has led to the computation of the orbits of seven new 
solar comets seen in the field-of-view (FOV) of the CP. The activities nec- 
essary to meet these large-angle attitude-pointing sequences, such as slew 
sequence planning, viewing-period prediction, and tracking-bias computa- 
tion are described. Analysis is presented for the computation of 
manuevers and pointing parameters relative to the SMM-unique, Sun- 
centered reference frame. Finally, science data and independent attitude 
solutions are used to evaluate the large-angle pointing performance. 
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1. INTRODUCTION 


This paper presents 

• Algorithms used in support of large-angle attitude pointing maneuvers by the 
Solar Maximum Mission (SMM) spacecraft to nonsolar scientific objectives 

• A brief history of the mission and its primary scientific objectives 

• Background of the problem of pointing at celestial objects, especially as pertains 
to the SMM 

The algorithms used to compute the necessary parameters for the observations are de- 
rived and results are presented that show the algorithms worked correctly, with accuracy 
significantly better than that specified by the SMM project. This work was performed by 
the Flight Dynamics Division (FDD) attitude determination and control analysts in the 
Flight Dynamics Facility (FDF) at Goddard Space Flight Center (GSFC). 

1.1 MISSION HISTORY 

The SMM spacecraft was launched in February 1980 from the Eastern Test Range into an 
approximately circular low-Earth orbit, with an inclination of nearly 28 degrees (deg) 
(Reference 1). The spacecraft functioned normally until November 1980 when the stand- 
ard reaction wheel (SRW) package that controls the spacecraft attitude failed. To preserve 
the mission, the spacecraft was put into a spin, approximately 1 deg per second (sec), 
about the minor principal axis. While the spacecraft was in this spin mode, very little 
scientific work was accomplished. In April 1984, the spacecraft was repaired in-orbit by 
astronauts aboard the Space Transportation System (STS) orbiter. The entire modular 
attitude control system (MACS) was replaced, and the spacecraft returned to normal 
operations. 

Some time before the repair mission, scientists at the GSFC Laboratory for Astronomy 
and Solar Physics proposed that the SMM be used to observe Comet Halley for approxi- 
mately 1 month before and after perihelion of the comet’s orbit (February 9, 1986). Dur- 
ing this time, the comet, as seen from the Earth, would be between 7 and 45 deg from the 
Sun, and very little useful ground-based observation was anticipated. 

1.2 SMM SCIENTIFIC OBJECTIVES AND PAYLOAD 

The original scientific objective of the SMM was to study solar phenomena, especially the 
solar maximum of 1981. The mission was planned to be three-axis stabilized to keep the 
scientific payload pointed at the Sun. However, because of the SRW failure, most of the 
scientific instruments missed the solar maximum. Since the repair of the spacecraft, a 
tremendous amount of solar radiation data at a variety of wavelengths has been amassed. 
In addition, beginning with the Comet Halley observations, several interesting nonsolar 
targets have been and are being studied. 

A Coronagraph/Polarimeter (CP), designed to study the solar corona in the visible spec- 
trum, was used to observe Comet Halley and the Moon. In addition, seven solar grazing 
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comets have been discovered during routine operation of the CP. A Hard X-Ray Burst 
Spectrometer (HXRBS) was pointed at Cygnus X-l (a suspected black hole) and the Crab 
Nebula. The Gamma Ray Spectrometer (GRS) was used to study gamma ray output from 
Supernova 1987A and the Ultra-Violet Spectrometer (UVSP) observed a transit of the Sun 
by Mercury. All of these observations required attitude excursions from the primary solar 
objective and provided valuable data for the scientific community on a wide range of 
celestial objectives. 

1.3 MANEUVER REQUIREMENTS 

The goal of this analysis is to provide algorithms for attitude control of the SMM so that 
the scientific payload points at the desired celestial body and tracks it, even if it is moving 
relatively fast, to keep the body in the instrument field of view (FOV). In addition, the 
periods during which the celestial body is visible to the spacecraft (i.e., not occulted by 
the Earth) must be determined. 

Requirements for attitude support of the large-angle attitude pointing maneuvers were not 
defined until several years after the spacecraft had been launched. For this reason, 
mathematical algorithms were developed for each nonsolar objective, as needed. The 
algorithms were later combined to form one software system that could efficiently plan 
the large-angle off-solar pointing maneuvers for any celestial body. The software was 
designed to plan the maneuver sequence, compute the viewing periods available, and 
determine the tracking parameters. Once computed, these data are uplinked to the space- 
craft. 

The constraints imposed by the SMM experimenters were few and were concerned with 
getting the best observations possible. To accomplish normal pointing, an accuracy of a 
tenth of a degree was desired. Another consideration was the determination accuracy in 
the viewing periods. Because of uncertainties in the spacecraft ephemeris, pointing accu- 
racy, and spacecraft gyro drift, accuracy in computing the start and stop times of the 
viewing periods was estimated to be 10 sec (3(7). 

Consideration for spacecraft safety placed other constraints on the maneuvers. The most 
important constraint was that the spacecraft had to remain in a power-positive mode to 
ensure survival of the spacecraft following any observation period contingency. This effec- 
tively meant that the body-fixed Solar Panel Array would always need to be absorbing at 
least as much energy as the spacecraft was using. During nominal operation, the solar 
panels are approximately perpendicular to the spacecraft-to-Sun vector. As a result of the 
analysis by the power engineers, the maximum off-Sun pointing angle was set at 65 deg. 

2. PROBLEM ANALYSIS 

2.1 THE SMM REFERENCE FRAMES 

The SMM spacecraft attitude is represented in a solar reference frame, referred to as the 
SUN frame. The X-axis of the SUN frame is defined by the apparent spacecraft-to-Sun 
unit vector, V ss . The Y-axis is the unitized cross-product of the Sun’s North Pole vector, 
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s„p , and . The Z-axis completes the right-handed coordinate system. The equations 
for the unit vectors along the Sun frame axes are (Reference 2) 


>1 

II 

XI 

(1) 

¥ = a (Snp x 20 

(2) 

a (Snp - (Sjip ' X) X) 

(3) 


where 

a = [1 - (&p • X) 2 r' /2 (4) 

Figure 1 illustrates the SUN coordinate system. 



Figure 1. Solar Coordinate System 


The transform, A, from geocentric inertial (GCI) coordinates to the SUN frame, is com- 
puted from the actual spacecraft and Sun positions with consideration for the Earth and 
spacecraft velocities and the speed of light to correct for aberration. Thus, for attitude 
determination, reference vectors must be transformed using A to get the reference vectors 
in the SUN frame. Due to the spacecraft’s orbit about the Earth and the Earth’s orbit 
about the Sun, the SUN frame is constantly changing. Thus, each time the attitude of the 
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spacecraft is computed by the spacecraft’s onboard computer (OBC) or the ground sys- 
tem, the SUN frame has to be recomputed. 

The SMM spacecraft attitude is represented as the transformation from the SUN frame to 
the SMM body frame. As shown in Figure 2, the fine-pointing Sun sensor (FPSS) and the 
scientific instruments are at one end of the spacecraft. During normal operations, the 
attitude control system is designed to point this end of the spacecraft at small areas of the 
Sun. The X-axis (roll) of the spacecraft body frame is defined to coincide with the direc- 
tion of the FPSS boresight. Two orthogonal axes perpendicular to the FPSS boresight and 
corresponding to the FPSS outputs are defined as the Y- (pitch) and Z- (yaw) axes of the 
body frame. This alignment is advantageous because the spacecraft pitch and yaw rota- 
tions are read directly from the output of the FPSS when the Sun is in the FOV of the 
FPSS, which is 4-square deg. Thus, during normal operations, the spacecraft is observing 
the Sun, and the attitude accuracy in pitch and yaw is 5 arc-sec (3(7). The roll is deter- 
mined by the fixed-head star trackers (FHSTs) to an accuracy of approximately 30 arc- 
sec (3(7). During periods of solar occultation or spacecraft pitch or yaw offpoint, the Sun 
is not in the FPSS FOV, and the attitude of all three axes is determined by the FHSTs to 
an accuracy of approximately 30 arc-sec (3(7). Most of the scientific instrument 
boresights are parallel to the boresight of the FPSS; thus, the output from the FPSS can 
be used directly to determine where the instruments are pointing to within 5 arc-sec when 
the Sun is being observed. 

The SMM attitude matrix, B, is the matrix that transforms a vector from the SUN frame 
to the SMM body frame. If the roll angle is 0, the pitch angle 6 , and the yaw angle f } , 
the matrix B is (Reference 2) 


B = T 3 (jS) TjOT T,«0 (5) 

where T* are the individual matrices for each Euler rotation on the intermediate axes. 
Thus, 

cos p cos 0 i cos P sin 0 sin 0 + sin p cos 0 1 - cos p sin 0 cos 0 + sin p sin 0 
B = - sin ft cos 0 1 - sin ft sin 0 sin 0 + cos P cos 0 | sin p sin 0 cos 0 + cos P sin 0 

sin 0 | - cos 6 sin 0 I cos 0 cos 0 

( 6 ) 

2.2 CALCULATION OF POINTING MANEUVERS 

Pointing at a celestial object involves rotating the SMM to align its X-axis with the 
spacecraft-to-target vector. In general, pointing from an initial target to any other target 
on the celestial sphere can be accomplished by consecutive rotations about any two 
orthogonal axes (Reference 3). Thus, any pair of rotations of roll, pitch, or yaw will 
suffice. For convenience, the first maneuver is chosen to be a roll. After this, the four 
choices for the second maneuver are positive or negative pitch or yaw. The magnitude of 
any of these second maneuvers is simply the angular distance from the Sun to the target 
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objective. For this algorithm, the roll corresponding to each choice of second rotations is 
computed, and the smallest roll is chosen. 


x, TO SUN 



Figure 2. Scientific and FPSS Boresight Directions 


To compute the maneuvers, the roll-yaw combination is computed first. The attitude ma- 
trix, B, for this case is T 3 (yS) T j(0) or 


B = 


cos ft 
-sin yS 
0 


sin yS cos 0 | sin yS sin 0 
cos yS cos 0 | cos /? sin 0 
- sin 0 | cos 0 


( 7 ) 


This transforms vectors from the SUN frame to the body frame after the maneuvers have 
been completed. 
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Since the scientific instrument boresights are parallel to the X-axis, the representation of 
the celestial target vector in the body frame after the maneuver is desired to be 


Vscb 



( 8 ) 


The GCI coordinate representation of the celestial target vector, unchanged by the maneu- 
ver is 


V 

JLgCI 


Vgcil 

Vgci2 

Vgci3 


( 9 ) 


To represent the celestial target vector in the SUN frame, the matrix A is calculated for 
the time of the maneuvers. The representation of the celestial target vector, Vs, in the 
Sun frame is 


Ys = A Vg C i 


( 10 ) 


while is related to V, as follows: 

V, = B T V, cb 


Vs, 


cos P | -sin ft | 0 


1 

V, 2 

= 

sin P cos <p | cos P cos <f> | - sin (j> 


0 

V s3 


sin P sin <p | cos P sin <p | cos 0 


0 


( 11 ) 


( 12 ) 


Thus, V* can be computed from Equation (10) and substituted into the left-hand side of 
Equation (12). The unknowns in Equation (12) are then the maneuver angles P and <p . 
Further reduction of this equation leads to 


V sl = cos p 

(13) 

V S 2 = sin p cos 0 

(14) 

V s3 = sin p sin <p 

(15) 
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Dividing Equation (15) by Equation (14) 


V S 3 /V s2 = tan 0 (16) 

Thus, the equations for the roll-yaw maneuvers are 

Yaw = p = cos' 1 (V sl ) (17) 

Roll = 0 = tan' 1 (V s3 /V s2 ) 08) 

The sign of the roll is resolved by using the function DATAN2 to compute the inverse 
tangent. The roll-pitch sequences are calculated from the above roll-yaw results by adding 
and subtracting 90 deg from the roll and selecting the sign of the pitch accordingly. As 
stated earlier, the maneuver sequence requiring the smallest roll rotation is the favored 
sequence. For targets that are further than 65 deg from the Sun, the maneuver to the 
attitude closest to the celestial body is computed. In this case the angular distance from 
the target to the planned closest attitude is also computed and delivered to the experi- 
menters. 

2.3 COMET POSITION DETERMINATION 

To date, seven Sun-grazing comets have been discovered in the CP FOV during normal 
observations of the solar corona. The CP uses an occulting disk on its sensor FOV to 
block out the main part of the Sun, enabling it to observe the solar corona. These comets 
could not have been observed without this special instrumentation of the CP. Several 
observations of each comet in CP instrument coordinates have been reported. In support 
of these observations, it is desired to transform the CP comet coordinate positions to GCI 
coordinates to determine the orbits of the comets. 

The determination of the comet positions in GCI coordinates is the inverse of the maneu- 
ver planning. It is assumed that the CP boresight is parallel to the X-axis of the space- 
craft. The comet position in CP coordinates is transformed to an equivalent spacecraft 
roll-pitch maneuver sequence that would point the spacecraft directly at the comet. The 
spacecraft, however, is never actually maneuvered; the maneuver is only computed so 
that the maneuver planning procedure can be reversed, and the comet position computed. 
The position vector in the Sun frame is then determined from the roll-pitch angles. The 
comet position vector in GCI coordinates is computed by multiplying the comet position 
vector in the Sun frame by the transpose of A at the time of the observation. 

The CP uses a two-dimensional polar coordinate system. The first coordinate is the coun- 
terclockwise rotation about the X-axis (roll) from the positive Z-axis (yaw). Thus, this is 
the negative roll necessary to align the positive Z-axis (yaw) with the Sun-to-comet vector. 
The second coordinate is the angular separation between the center of the Sun and the 
comet, along the positive Z-axis (yaw). This corresponds to a negative pitch angle. The 
result is the roll-pitch angle pair that would be required to maneuver the spacecraft and 
align the boresight with the spacecraft-to-comet vector. 
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If the matrix, B\ for a roll (0)-pitch ( 9 ) maneuver from a null attitude in the Sun frame 
to the comet is 


cos 9 | 

sin 9 sin 0 | 

- sin 9 cos 0 

0 1 

COS 0 | 

sin 0 

sin 9 | 

-cos 9 sin 0 | 

cos 9 COS 0 


(19) 


then the position vector, V c , of the comet in the Sun frame is simply 


or 


Yc 


= B 


'T 


1 

0 

0 



cos 9 

0 1 

sin 9 


1 

Yc = 

sin 9 sin 0 

COS 0 I 

-cos 9 sin 0 


0 


-sin 9 cos 0 | 

sin 0 | 

cos 9 cos 0 


0 


( 20 ) 


( 21 ) 


This assumes that the spacecraft attitude is zero roll, pitch, and yaw so that the SUN 
frame is aligned to the body frame. If the spacecraft is at a non-nominal attitude, the 
matrix B must be recomputed, taking into account the spacecraft attitude at the time of 
the observation. The comet position in GCI coordinates, V cgC j, is then 


V . - A t V 

JLcgci — ^ JLc 


( 22 ) 


The vector, V cgci , is then reduced to right ascension and declination (mean of 1950) and 
delivered to the International Astronomical Union (IAU), where it is used for orbit com- 
putations. 

2.4 CELESTIAL BODY VIEWING PERIOD PREDICTION 

Because of the geometry of the SMM orbit, every point in the celestial sphere of the 
spacecraft is occulted by the Earth at some time in one revolution of the ascending node 
of the spacecraft’s orbit. Thus, a list of periods when the celestial target can be viewed by 
the spacecraft is necessary for all times that observations of nonsolar targets are being 
planned. The computation of these periods is done by comparing the relevant angles in 
the geometry of the spacecraft, Earth, and the target objective. 

The unit vector, , from the center of the Earth to the spacecraft, is constantly chang- 
ing as the spacecraft orbits the Earth and is known from the spacecraft ephemeris. The 
vector from the Earth to the target body, Web. is normally computed from the right 
ascension and declination of the target body in the Earth’s celestial sphere. To compute 
the viewing periods, the spacecraft-to-target vector, W < h is needed. Since the distance to 
the viewing target in most cases has not been relevant, the Earth-to-target vector is 
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represented as a unit vector. Thus, vector addition of the Earth-to-target vector and the 
Earth-to-spacecraft vector does not give the spacecraft-to-target vector. Instead, it is as- 
sumed that the parallax of the target body’s position due to the spacecraft’s orbit around 
the Earth is negligible, and, therefore, 


W eb - W, b (23) 

For objects out of the solar system and the outer planets, this approximation will yield 
results that are well within the constraints. For the inner planets and comets, this approxi- 
mation leads to an error in the vector of between 1 arc-sec and 20 arc-sec. This will 
translate to an error of less than 1 sec in the occultation times, which can be considered 
negligible compared to the 10-sec constraint on the accuracy. If observations of the Moon 
are planned, the position error caused by the above assumption can reach 1 deg. This will 
translate to approximately a 15-sec error in the occultation times, which is greater than 
the maximum error allowed. To resolve this, the Earth-to-Moon vector obtained from 
lunar ephemerides is used to calculate W, b directly by vector addition. 

Once Ws b and Ws Ce are known, the angle, <£, between the vectors is computed by taking 
the inverse cosine of the dot product of the two vectors as 


$ = cos’ 1 (Wsb • Wsce) 


(24) 


This geometry is illustrated in Figure 3. 


SPACECRAFT 



Figure 3. Occultation Geometry 
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To compute whether the celestial target is blocked by the Earth, <t> is compared to the 
subtended Earth half angle, Q. The computation of Q assumes a perfectly circular 



Earth, with a radius equal to the equatorial radius. This translates to a maximum of 8 sec 
for celestial objects with declinations near positive or negative 90 deg. However, nearly 
all of the celestial targets observed by the SMM have had declinations close to the Sun s. 
Thus, the typical error induced by the Earth’s oblateness is typically near 3 or 4 sec. 

The half angle subtended by the Earth is given by 

Q = sin -1 (R/D) (25) 

where R is the equatorial radius of the Earth and D is the magnitude of the Earth-to- 
spacecraft vector obtained from the spacecraft ephemeris file. The relationship of the 
computed angles is also shown in Figure 3. The only variable on the right-hand side of 
Equation (25) is D. This varies with the orbital eccentricity and more slowly with the rate 
of decay of the spacecraft’s altitude above the Earth. Thus, for SMM, Q changes very 
slowly and is approximately 69 deg. The celestial object is considered to be occulted by 
the Earth from the spacecraft if Q is larger than 4>. Conversely, if T> is larger than Q , 
the celestial object is considered in view of the spacecraft. This makes sense intuitively, 
since at times the spacecraft is between the Earth and the celestial object, will be 
obtuse, while if the Earth is between the spacecraft and the celestial object, O will be 
acute. 

Earth occultation computations are performed every 4 minutes (min) for the period of 
interest. When a change occurs, a binary search is used to narrow the change time to the 
nearest second, which will typically have a maximum error (3 o') of less than 10 sec 
(Reference 4). 

2.5 TRACKING PARAMETER COMPUTATION 

Targets not in the solar system, such as distant stars, do not move significantly relative to 
the SMM due to the spacecraft’s orbit about the Earth or the Earth’s motion about the 
Sun. Thus, during observations of very distant objects, provisions need not be made to 
track motion of the celestial target due to parallax and proper motion. However, bodies in 
heliocentric (planets or comets) or geocentric (the Moon) orbits will move through an 
angle d , which can be represented as a rotation about a unit vector, e, in the spacecraft 
inertial reference unit (IRU) coordinate frame as illustrated in Figure 4. This motion is 
due to the spacecraft’s motion about the Earth, the Earth’s motion about the Sun, and the 
target objects motion in its orbit about the Sun or Earth. During an observation period at 
a constant attitude, this motion will cause the celestial target to move through, and even- 
tually out of, the scientific sensor FOV. To compensate, the spacecraft is rotated about 
the axis e by the angle 5 from the beginning to the end of the observation period. This 
motion is called tracking. 

Tracking by the spacecraft is accomplished by computing the angle, 5 , and axis of rota- 
tion, e, of the target body in the IRU sensor frame of reference and combining <5 and e 
with the observation time At to form a command gyro rate bias vector (historically called 
BGDT) which is uplinked to the spacecraft. These BGDTs cause the gyros to show motion 
opposite to the way the target is moving. This apparent (in that the spacecraft is not 
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actually moving) motion causes the SRWs to counteract the apparent motion and move 
the spacecraft in the direction of the celestial objective. Thus, tracking is accomplished. 
To compute the BGDTs, the transformation matrix, H, between the SMM body-frame 
attitude at the start of the observation period and the attitude at the end of the observation 
period is computed and used to determine average rates about the three IRU axes. 


TARGET AT 
OBSERVATION 
STOP TIME 



Figure 4. Tracking Geometry 


To compute H, the maneuvers to point the spacecraft at the target are computed by the 
methods described earlier at the observation start and stop times. The maneuvers are 
assumed to be relative to an attitude of roll, pitch, and yaw equal to zero. Thus, using 
Equation ( 6 ), the SUN-frame-to-body-frame matrices, Bi and B2, at the start and stop 
times are computed. However, the coordinate frame that B is referenced to at all times on 
the SMM mission is the SUN frame, which, as stated above, is not inertial. The attitude 
motion for target tracking must be referenced to an inertial frame. For this, the GCl-to- 
SUN frame matrix, A, is used. In addition, because the BGDTs are for the IRUs, they 
must be represented in the IRU sensor frame. Defining Gi and G2 to be the GCI-to-IRU 
transformation matrices at the start and stop times, respectively, and M to be the IRU-to- 
body transformation matrix, then 

Gj = M t Bj A] (26) 
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G 2 = M t B2A 2 


(27) 


and 


G 2 = HG, 


(28) 


Solving for H yields 


H = G 2 GT = M T B 2 A 2 (M t B, A0 t 


= M t B 2 A 2 Aj Bj M 


(29) 


If the Sun frame were an inertial frame, Ai would be equal to A 2 , and Equation (29) 
would be simplified. Now that H has been found, <5 and £ can be computed from H, and 
the BGDTs can be computed. 


To compute d and e, H is written as an Euler rotation with unit vector e in the IRU frame 
and rotation angle <5 . Thus, using the Euler convention for a general rotation, H is (Refer- 
ence 4) 


H = cos ((5) I + (1 - cos (<5)) ee T - sin (d) [[e]] (30) 


where [[e]] is the skew symmetric matrix representation of a vector. The general form for 
a three-component vector is 


[lx]] 


0 

-x 3 

x 2 

x 3 

0 

-Xi 

-X 2 

Xi 

0 


(31) 


Thus, it is desired to manipulate Equation (30) to solve for e and <5 . Defining the pa- 
rameter k 2 as 


k 2 = (H 23 - H 32 ) 2 ♦ (Hj, - H,j) 2 + (H 12 - H 2 ,) 2 ( 32 ) 


substitution from Equation (30) leads to 

k 2 = [e 2 e 3 (1 - cos 6) + e 3 sin <5 - e 2 e 3 (1 - cos 6) + ti sin d} 2 


+ [ei e 3 (1 - cos d) + e 2 sin <5 - e 3 e 3 (1 - cos (5) + e 2 sin (5] 2 (33) 
+ [e 3 e 2 (1 - cos 6) + e 3 sin <5 - e 3 e 2 (1 - cos 5) + e 3 sin d] 2 
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Some simple algebra leads to 


k 2 = At\ sin 2 <5 + 4e^ sin 2 6 + 4 sin 2 6 


Since e is a unit vector, the 
to 1. Therefore, 


(34) 

4 * sin 2 d * (e 2 + t\ + e 2 ) 

last quantity on the right-hand side of Equation (34) is equal 


k 2 = 4 * sin 2 6 

(35) 

or taking the square root of both sides and solving for 6 


d = sin ' 1 (k/ 2 ) 

(36) 

where k is defined by Equation (32). 


To solve for £, the off-diagonal terms in Equation (30) are used, 
using ei. 

To demonstrate this 

H 23 - H 32 = e 2 e 3 (1 - cos 6) + ei sin 6 - e 2 e 3 (1 - cos 

6) + e, sin 6 

= 2 * ej * sin 6 

(37) 

Solving for ei 


e i = (H 23 - H 32 )/(2 * sin <5) 

(38) 

Similarily, 


e 2 = (H 31 - H j3 )/(2 * sin 6) 

(39) 

e 3 = (H ,2 - H 2 i)/( 2 * sin <5) 

(40) 

To compute the BGDTs, the angle of rotation in the IRU coordinate 
the observation time. The rate about each of the IRU axes is then 

frame is divided by 

x' = ej * ((5/At) 

(41) 

y' = e 2 * ((5/At) 

(42) 

z' = e 3 * (6/ At) 

(43) 


where At is the observation time. These rates are converted to the units of the SMM IRU, 
radians per 0.512 sec, and uplinked to the spacecraft for use in the control loop. 
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3. RESULTS AND CONCLUSIONS 


3.1 RESULTS 

To date, these algorithms have been used to support all the SMM scientific off-Sun ma- 
neuvers. Since the computed parameters are designed to support pointing to celestial 
objects other than the Sun, the performance of the algorithms is ultimately measured by 
the quality of these observations. Thus far, the science observations have been completed 
without incident and have yielded excellent results each time. Images obtained from the 
SMM scientific instruments have always had the target object in the sensor FOV and at 
the correct times. In addition, for the heliocentric and geocentric objectives, the tracking 
system has continuously kept the object in view. 

As an independent check, attitude determination solutions during the observations have 
indicated the spacecraft pointing accuracy is within 180 arc-sec. The accuracy of the solu- 
tions was approximately 30 arc-sec (3 cr) . The bulk of the pointing error was due to uncer- 
tainty in the IRU bias. During the observations, the spacecraft was inertially stabilized by 
the IRU package. However, spacecraft attitude drift occurred because of uncertainty in 
the knowledge of the IRU bias. This bias was calibrated periodically, and the error in the 
pointing did not exceed 120 arc-sec due to the error in the bias. 

Most of the remaining error in the pointing accuracy was due to uncontrollable timing 
problems with the maneuvers. Since the geometry of the spacecraft, Earth, Sun, and 
target-body system was constantly changing, the required maneuvers were constantly 
changing. The algorithm derived above assumes the maneuvers were instantaneous. How- 
ever, the SMM maneuvered at approximately 5 deg per min; consequently, pointing 
errors were incurred. In addition, the start time of the maneuvers was calculated at 5-sec 
intervals. Thus, if the maneuvers did not start at these times, greater error was incurred. 
The pointing errors due to the timing errors did not exceed 50 arc-sec. Other potential 
error sources include uncertainties in the position of the target body, the initial spacecraft 
attitude from which the maneuvers were performed, and the IRU scale factors which 
controlled the magnitude of a commanded maneuver. 

3.2 CONCLUSIONS 

Algorithms to support off-Sun maneuvers of the SMM have been presented. These algo- 
rithms compute parameters for maneuver planning, viewing period start and stop times, 
and target tracking rate. All of these items are necessary to maneuver the spacecraft to 
observe and track the celestial target objective. Because the Sun reference frame is non- 
inertial, the equations differ from ones normally seen. This work can be adapted to space- 
craft that are controlled relative to the inertial frame; however, the primary application is 
to missions with non-inertial reference frames, such as Earth-pointing (1 revolution per 
orbit) missions. 
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